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Revisiting the Impact of Axions in the CooUng of White Dwarfs 
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Abstract. It has been shown that the shape of the luminosity function of white 
dwarfs can be a powerful tool to check for the possible existence of DFSZ-axions. In 
particular, Isern et al. (2008) showed that, if the axion mass is of the order of a few meV, 
then the white dwarf luminosity function is sensitive enough to detect their existence. 
For axion masses of about nia > 5 meV the axion emission can be a primary cooling 
mechanism for the white dwarf and the feedback of the axion emission into the thermal 
structure of the white dwarf needs to be considered. Here we present computations of 
white dwarf cooling sequences that take into account the effect of axion emission in a 
self consistent way by means of full stellar evolution computations. Then, we study and 
discuss the impact of the axion emission in the white dwarf luminosity function. 



1. Introduction 



X: 



The Peccei-Quinn mechanism is one of the most convincing explanati ons for the ab - 
sence of CP- violating effects arising from the QCD vacuum structure (see lRafFellll996h . 
A natural conse quence of this mechanism is the existence of a new particl e, a boson 
called the axion ( Wilczekiri978h . One of the axion models, the DFSZ model dZhitnitsky 
19801 : iDine et al 11981 ) allows for the interaction of axions with charged leptons. The 
existence of s uch DFSZ- axio ns would have an impact in the cooling of white dwarfs as 
pointed by Raf feltl(ll986h and llsem et all ( 19921 ) . Because the evolution of white dwarfs 
is mostly a simple cooling process and the basic physical ingredients needed to predict 
their evolution are relatively well known, white dwarfs offer a unique oportunity to 
test new physics under conditions that can not be obtained in present day laboratories. 
Isern et aU (2008) s howed that \ yith the current knowled ge of the white dwarf luminos- 
ity function (H arris et al.ll2006l : iDe Gennaro et al.ll2008 ) it might be possible to detect 
axions as light as ~ 5meV. For axion masses which might be detectable through the 
white dwarf luminosity function (m^ > 5meV) t he axion cooling is comparable to the 



neutrino and photon cooling of the white dwarf dlsem et alj|200 8'). In such a situation, 
we expect the existence of a signifficant axion emission to impact the thermal structure 
of the white dwarf, and consequently to alter both the photon and neutrino emission. 
Thus, for the range of interest of the axion masses we expect that the axion emission 
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Figure 1. Left: Comparis on of the axion {n tq = 5meV) emission of our 0.609Mq 
and the 0.61Mq sequence ofl lsern et alj ( 120081) . The effects of the departure from the 
isothermal core approximation can be apreciated at Mboi < 6. Right: Axion (black 
curves) and neutrino (blue curves) emission for our O.6O9M0 sequences for different 
axion masses. The impact of the axion emission in the thermal structure of the white 
dwarf can be appreciated in the decrease of the neutrino emission at higher axion 
masses. Clearly, axion emission can not be treated perturbatively at ntg > 5meV. 



can not be treated as a perturbation to the white dwarf cooling and that a self consistent 
treatment of the axion emission is necessary. In the present work, we improve previous 
works by studying the impact of the axion emission in the cooling of white dwarfs by 
means of a self consistent treatment of the axion emission and state of the art white 
dwarf models. 



2. White dwarf models, input physics and the white dwarf luminosity function 



The initial white dwarf models were taken from Renedo et al. ] i20\(f). Specifically 4 



different initial white dwarf models of 0.524Mq, 0.609Mo, 0.705Mq and 0.877 Mq 
were adopted. These models were obtained from computing the complete evolution of 
initially IMq, 2Mq, 3Mq and 5 Mq ZAMS stars with Z-O.Ol, which is agreement with 
semi-empirical determinations of the initial-final mass relationship ( Salaris et al.,.2009.) . 
It is worth noting that these sequences were computed from the ZAMS to the thermal 
pulses at the AGB and finally to the post-AGB and white dwarf stages. For each ini- 
tial white dwarf model 6 cooling sequences with different assumed axion masses were 
computed {ma= 0, 5, 10, 15, 20 & 30 meV). Axion emission by both Compton and 
Bremsstrahlung processes was included, although in white dwarfs only the latter is 
important. Axion Bremsstrahlung emission under degenerate conditions was included 
adopting the prescriptions of Nakaaawa et al. ( 1987, 1988) for strongly coupled plasma 
regime (T > 1) and iRaffelt & We iss ( 1995) for weakly coupled plasmas (F < 1). Com- 
putations were performed with LPCODE stellar evolution code, which is specifically tai- 
lored for the computation of white dwarfs and includes all relevant microphysics such 
as equation of state, radiative and conductive opacities, element diffusion and can even 
handle the effects of phase separation, crystallization and release of latent heat. 
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Figure 2. White Dwarf luminosity functions constructe d for the diff erent axion 
masses compared with the luminosity functions derived by [Harris et alj ()2006) and 
iDe Gennaro et alj (|^08). DFSZ-axions heavier than > lOmeV are clearly exp- 
cluded by the observed white dwarf luminosity functions. 



T o construct theoretical white dwarf lu minosity functions we f ollowed llsem et al. 
( 2008h but adopted the method presented by llben & LaughlinI ( 19891) . From the cooling 
times, tcil, m), computed with LPCODE for each value of the axion mass, we computed 
the number density of stars per luminosity as 



(1) 



For a given white dwarf luminosity bin (Z), and mass of the progenitor (M), the forma- 
tion time of the star, t, is given hy t + tsviM) + tdl, m) = T^. To compute the white dwarf 
luminosity function we adopt the following additional ingredients: A Salpete r initial 
mass function N{M), the initial-final mass relatioi iship m(M) from.Salaris et al.l (12009.) . 
the stellar lifetimes fev from the BaSTI database (iPietrinfemi et al. 2004) and constant 
star format i on rat e ip. Normalization of theoretical luminosity functions is done as in 
Isemetal.l(l2008h . 



2.1. Results and future work 

In Fig. [T] (left panel) the axion emi ssion of our 0.609Mo sequence is compared with 
the O.6IM0 sequence of llsem et al.l (i2008.) for an axion of nia = 5meV. There is an 
overall good agreement between both predictions at low luminosities (Mboi > 6). The 
departure between both curves at high lumino sities (Mboi < 6) can be traced back to the 
isothermal core approximation of llsem et all (1200 8) which leads to an underestimation 
of the axion emission at high luminosities when the maximum temperature of the core 
is located off-centered. The right panel of Fig. [T] shows both the axion and neutrino 
emissions of our 0.609Mo sequences computed with different axion masses. As ex- 
pected the higher the axion mass, the higher the axion emission and the cooling speed. 
A more interesting feature can be seen by analysing the neutrino emission at different 
axion masses. As can be seen in Fig. [T]the neutrino emission is reduced as the axion 
mass is increased. This result is due to the fact that when axion emission is included 
this leads to an extra cooling of the white dwarf core which alters the thermal structure 
of the white dwarf (as compared with the case with no axion emission) which, in turn. 
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leads to a decrease of the neutrino emission at a given surface luminosity of the star. 
Note that this is true even for our lighter computed axion mass (nia = 5meV), for which 
the axion emission is already different from the case with no axions (nia = OmeV). 
Then, the inclusion of the feedback effects of the axion emission on the thermal struc- 
ture of the white dwarf, leads to a decrease of the neutrino emission which will dimin- 
ish the sensitivity of the white dwarf cooling times to the axion mass. Consequently 
the axion emission needs to be treated self-consistently when dealing with axions in 
the range detectable through the white dwarf luminosity function. In Fig. |2]we show 
the resulting white dwarf luminosity functions for each axion mass as compared with 
the observationally derived ones. It can be clearly seen that axion masses larger than 
10 meV would lead to strong disagre ements wit h the luminosity functions derived by 
Harris et al. (2006) and De G ennai o et al.l (l2008h . On the other hand, the existence or 
not of DFSZ-axions with masses lower than 5 meV can not be concluded without a 
detailed statistical analysis of the uncertainties. 

We have shown that the impact of the axion emission on the white dwarf lumi- 
nosity function will be overestimated if the effect of axion emission on the thermal 
structure of the white dwarf is not taken into account. This is particularly important at 
log Ly > - 1 . Our preliminary analysis shows that the observed luminosity functions 
are consistent with the absence of any additional cooling mechanism (like axions), al- 
though a detailed statistical analysis of the results should be made before making a 
final statement. On the other hand, DFSZ-axion masses larger than nia > 10 meV, are 
clearly excluded by our present knowledge of the white dwarf luminosity function, in 
agreement with llsem et all (l2008h . 
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